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ABSTRACT 

The luminous Young Stellar Object (YSO) IRAS 07422-2001 is studied in the infrared. We discover star forming activity in embedded 
clusters located in a cloud detected at mid-IR wavelengths in emission. Multiple outflows are discovered from these clusters in the H2 
ro-vibrational line at 2.122yL/m. We detect at least six outflows from the cluster associated with the IRAS source and another outflow 
from a source located in a cluster detected ~2.7 arcmin NE of the IRAS source. Additional star formation is taking place in two other 
cluster candidates within the cloud. Three of the YSOs in the cluster associated with the IRAS source are detected at W.lfim at an 
angular resolution of ~0.8''. We have a tentative detection of a circumstellar disk in this cluster, seen as an extinction lane in the / 
and //-band images. The spectral energy distributions (SEDs) of the dominant YSOs in the cluster associated with the IRAS source 
and in the NE cluster are studied using radiative transfer models and the properties of the YSOs are estimated. The YSO associated 
with the IRAS source is probably in a very early Class I stage of formation. The source identified as the dominant YSO in the NE 
cluster appears to be older than the dominant YSO in the cluster associated with the IRAS source, but its observed flux seems to 
be contaminated by extra emission, which suggests the presence of a young source contributing to the SED at far-IR wavelengths. 
The star formation observed in the field of IRAS 07422-2001 supports the idea of hierarchical formation of massive star clusters and 
the growth of massive young stellar objects near the centres of multiple sub-clusters in a star forming clump through competitive 
accretion. 
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1. Introduction 

Studies of massive star formation have been getting significant 
attention in the recent years. The intriguing problem is to un- 
derstand how massive stars form - whether it is through a pro- 
cess similar to that for low mass stars, ie., through disk accre- 
tion and driving molecular outflows (e.g. Yorke & Sonnhalter 
120021) or through alternative scenarios like coalescence of lower 
mass stars (e.g. Bonnell, Bate & Zinnecker 1998 ). Even though 
it is still debatable how very massive stars form, it is almost 
clear that stars at least up to late-0 spectral types form primarily 
through disk accretion and that the poor collimation of outflows 
from many of these objects previously seen in single-dish CO 
observations is main ly due to mult iple co llimated outflows (eg. 
Beuther et al. 120021 Varricatt et al. l2QlQI) . However, it has been 
proposed that their formation may be diflTerent from low mass 
stars' in the sense that they may form mainly in massive clusters 
through competitive accretion, where the cluster potential aids 
the growth of massi ve YSO s near the centres of massive clusters 
(eg. Bonnell & Bate [2QQ6l) 

IRAS 07422-2001 (henceforth IRAS 07422) is a lumi- 
nous YSO, first detected in i2cO(J=l-0) by Wouterloot & 
Brand (1989 ) using the 15-m SEST. They detected a blue wing 
emission in CO, implying the presence of outflows. From a 
V/^^=52.7km s"^ of the CO line, they derived a kinematic dis- 
tance of 5.1 kpc and a luminosity of 7.4xlO^L0. Palla et al. 
( 1993 ) detected strong H2O masers at ypeak=^^-^ km s"^ and de- 
rived a kinematic distance of 4.2 kpc. Brand et al. (119941 ) also 
detected H2O masers from this source at similar radial velocity 
and flux density. 

Slysh et al. (119991) did not find any 6.7 GHz class-II methanol 
maser emission from this region at a 3cr sensitivity of 3 Jy sug- 



gesting that the YSOs here may be very young. Codella & Felli 
(1995 ) did not detect any H89ck hydrogen recombination line 
emission from IRAS 07422 suggesting that this source is very 
young and has not yet formed a diflTuse Hii region. Urquhart et 
al. (1^507 1). in their Red MSX Source (RMS) radio survey, did 
not detect any emission (above their rms noise levels of 0. 14 mJy 
and 0. 17 mJy at 3.6 cm and 6 cm respectively) from IRAS 07422 
indicating that the YSOs detected here at mid- and far-IR wave- 
lengths are in pre-UCHii phase. 

A dense core was detected towards this source in 97.981- 
GHz CS(2-1) emission (Bronfman, Nyman & May 1996) at 
V/^r=52.8 km s"^ Slysh et al. ( 1997 ) observed OH thermal emis- 
sion at a similar velocity (52.3 km s"^). These are significantly 
diflTerent from the velocity at which H2O masers were detected. 
However, H2O masers seen towards massive YSOs are under- 
stood to be excited by jets (eg. Felli, Palagi, Tofani 1992, Goddi 
et al. 2005 ) . Hence we use the 52.8 km s"^ of the CS emission for 
estimating the distance. Using the latest rotation curves of Reid 
et al. (I2QQ9L we estimate a kinematic distance of 4.01 (+0.69, - 
0.64) kpc towards IRAS 07422. We will use this distance for the 
calculations in this paper. 

Even though outflow signatures were observed in the data 
of Wouterloot & Brand (119891) , no systematic observations have 
been reported on this source to investigate it further. Near- 
IR imaging is a powerful tool to study star forming regions. 
Outflows from YSOs emit strongly in the near-IR ro-vibrational 
lines of H2; the H2 v=l-0 S(l) line at 2.1218yum is a powerful 
tracer of outflows via shock excitation. The high spatial reso- 
lution available with moderately sized telescopes at IR wave- 
lengths when compared to that available with single-dish radio 
and mm telescopes and the much lower interstellar extinction at 
IR wavelengths compared to that at optical wavelengths make 



1 



Watson R Varricatt: Multiple outflows from the luminous YSO IRAS 07422-2001 



Table 1. Log of WFCAM observations 



UTDate 


Filters 


Exp. time 


Total int. 


FWHM 


(yyyymmdd) 


used 


(s) 


time (s) 


(arcsec) 


20110218 


IH, 


10,5, 


720, 360, 


0.83,0.81, 




K,H2 


5,40 


360, 800 


0.82, 0.81 


20110309 


K,H2 


5,40 


360, 800 


1.40, 1.28 


20110310 


H2 


40 


800 


1.16 


20110311 


H2 


40 


800 


1.37 



this regime ideal for studying massive star forming regions. In 
this paper, we use near-IR imaging to understand the star forma- 
tion in IRAS 07422. The new data are used along with available 
archival data to study the properties of the YSOs. 

2. Observations and data reduction 

2.1. WFCAM imaging 

Observations were obtained with the 3.8-m United Kingdom 
Infrared Telescope (UKIRT) and the Wide Field Camera 
(WFCAIVQ Casali et al. 2007 ) using near-IR broadband /, H 
and K MKO filters and a narrow-band MKO filter centred at 
the wavelength of the H2 (1-0) SI line at 2.1218yum. WFCAM 
has a pixel scale of 0.4'' pix"^ and employs four 2048x2048 
HgCdTe HAWAII-2RG arrays, each with with a field of view of 
1 3.65' X 13.65'. The four arrays are arranged in a square pattern 
in the focal plane, with a gap of 12.83' each in the fields cov- 
ered by adjacent arrays. For the current work, we used the data 
from only one of the four arrays in which the object was located. 
Observations were performed by dithering the object (to 9 points 
in the focal plane for /, H and K bands and to 5 points for H2) 
separated by a few arcseconds. We used a 2x2 microstepping, 
hence the final pixel scale is 0.2" pixel'^ The spatial resolution 
is limited by the seeing at the time of observations. 

The first set of observations were obtained on 201 10218 UT. 
Since outflows were detected in the H2 line, we further observed 
the region in H2 and K filters on three more epochs to improve 
the depth of detection. Table [T] gives a log of the WFCAM ob- 
servations performed. 

The data were reduced by the Cambridge Astronomical 
Survey Unit (CASU); the archival and distribution of the data 
are carried out by the Wide Field Astronomy Unit (WFAU). The 
sky conditions were not photometric during the observations. 
However, since the flux calibration was performed by deriving 
the photometry of isolated point sources with 2MASS detection 
in the large field of each array, the photometric accuracy is good 
for the stars detected. The derived magnitudes are in the MKO 
photometric system. Fig. [T] shows a JHH2 colour composite im- 
age of a 5.5'x5.5' section of the full WFCAM image, showing 
the area where the outflow and star formation activity discussed 
in this paper take place. 

The sky background was fitted and removed from the mo- 
saics. The H2 images were continuum- subtracted using the K- 
band images. The average of the ratio of counts, for a few iso- 
lated point sources, between the K and H2 images was obtained. 
The background- subtracted ^-band image was scaled by this ra- 
tio and was then subtracted from the background- subtracted H2 
image. The H2 mosaic of 20110218 was continuum- subtracted 
using an average of the two ^-band mosaics obtained on the 

^ http://www.j ach.hawaii.edu/UKIRT/instruments/wfcam/ 



same night; those of 20110309-20110311 were continuum- 
subtracted using the average of the two ^-band mosaics ob- 
tained on 20110309. Due to diff'erences in seeing between the 
K and H2 images, point sources give positive and negative resid- 
uals. A uniform continuum- subtraction over all objects in the 
field is difl&cult to achieve in regions with a wide range in near- 
IR colours. The K/112 ratio has a dependence on the apparent 
near-IR colours; hence highly reddened sources and objects with 
excess appear with large negative residuals in the continuum- 
subtracted H2 image (see Varricatt .20 11 i f or a discussion). Fig.|2] 
shows the continuum- subtracted H2 image of the region shown 
in Fig. [T] Images from all four epochs of our H2 observations 
are averaged here to improve the S/N of the emission features 
detected. 



2.2. Archival data 

The IRAS mission detected a point source at 0^=7:44:27.9 S=- 
20:08:32 (J2000) in afl four bands (at 12, 25, 60 and 100 yum). 
IRAS 07422 was detected well in the all sky survey conducted 
by the AKARI satellite (Murakami et al., [20071). The mid-IR 
camera - Infrared Camera (IRC) - on board AKARI detected 
three sources at 1 8 jim. The brightest two of these were also de- 
tected by the AKARI Far-Infrared Surveyor (FIS) in all 4 bands 
(at 65, 90, 140 and 160 yum). 

At mid-IR wavelengths, the MSX mission (in four bands 
A, C, D and E centred at 8.28, 12.13, 14.65 and 21.23 yum re- 
spectively) also detected three sources in this region. The loca- 
tions of the MSX sources roughly coincide with those of the 
three sources detected by AKARI-IRC at 18yum. This region 
was observed in the Wide-field Infrared Survey Explorer (WISE; 
Wright et al. l2QT0j mission in four bands - Wl (3.4yum), W2 
(4.6yum), W3 (12/jm) and W4 (22yum). The images in the four 
bands have angular resolution of 6.1, 6.4, 6.5 and 12.0" respec- 
tively and have an astrometric precision better than 0.15"for 
sources detected with good S/N. The excellent sensitivity of the 
WISE data enables detection of many faint features in the data. 

Most of the mid-IR observations available on this source 
have poor angular resolution compared to our near-IR WFCAM 
images. No Spitzer data is available for IRAS 07422. The only 
mid-IR observation available with good angular resolution is 
with UKIRT and Michelle. Michelle (Glasse et al. [T9971) is 
a mid-infrared imager/spectrometer operating in the 8-25 yum 
wavelength regime, using an SBRC Si:As 320x240-pixel ar- 
ray. It has an image scale of 0.21"/pix and a field of view 
of 67.2"x50.5"at the Cassegrain focal plain of UKIRT. The 
Michelle observations were performed on 20040422 (UT) using 
an filter centred at 1 1.2 yum 

The Michelle data were downloaded from UKIRT data 
archive and were reduced using the facility data reduction 
pipeline ORACDR. The observations were performed by chop- 
ping the secondary and nodding the telescope by 20" each (peak- 
peak) in mutually perpendicular directions to detect faint sources 
in the presence of strong background emission at thermal wave- 
lengths. The resulting mosaic has four images of the source (2 
positive and 2 negative beams). The final image was constructed 
after negating the negative beams and combining all 4 beams. 
Fig. [8] shows the final mosaic of the region around IRAS 07422. 
An average FWHM of 0.8" was measured from the sources de- 
tected in the 11.2 yum image. 
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Fig. 1. WFCAM JHH2 colour-composite image (/-blue, //-green, //2-red) in a 5.5'x5.5' field covering the embedded clusters 
associated with IRAS 17527. 



3. Results and discussion 

Our near-IR images reveal a very active picture of star formation 
associated with IRAS 07422. Fig. [T] shows a large number of 
red sources in at least two embedded clusters and two additional 
cluster candidates active in star formation. Regions containing 
them are enclosed in dashed boxes or ellipses in the figure and 
are labelled "CI-CIV"; the most prominent of these are "CI" 
and "CII". Interesting point sources are labelled "1-11" on the 
figure. 

The continuum- subtracted H2 image (Fig. |2]) shows mul- 
tiple outflows from this region. Dashed arrows are shown on 
Fig. [2] where aligned H2 line emission knots are present, which 



suggest the directions of the outflows. The H2 emission fea- 
tures are named MHO 1412-1420 in the Catalogue of Molecular 
Hydrogen emission-line Object^ (Davis et al. 2010 ) and are la- 
belled on the figure. Table |3] shows the directions and lengths of 
the outflows detected in Fig. [2] The angles given in the table are 
east of north, measured from our H2 image. 

Near-IR (J - H) - (H - K) colour-colour diagram is a pow- 
erful tool to identify YSO candidates via reddening and ex- 
cess. The JHK magnitudes of point sources detected over the 
13.5'xl3.5'field covered by the array on which IRAS 07422 
was located are used to construct the colour-colour diagram. 



2 http://www.jach.hawaii.eduAJKIRT/MHCat/ 
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Fig. 2. The continuum-subtracted 2.122-yum image of IRAS 07422 over a 5.5'x5.5' field shown in Fig.[l] smoothed with a 2-pixel 
FWHM Gaussian to enhance the appearance of the faint emission features, "x" and show the locations of the IRAS and MSX 
sources respectively. The top-right inset shows an expanded view of a TxT field containing "CI"; the contours generated from the 
11.2-yum Michelle, showing the 3 objects detected at 1 1.2yum (#1, #2 and #3), are overlaid on it. The top-left inset shows a 35''x35'' 
field in "CII". 



Magnitudes of the sources detected over multiple epochs are 
averaged. We had clouds present during the observations. 
Therefore, to avoid any influence from variable clouds, objects 
located within 7'' of the edges of our near-IR images are omit- 
ted from the colour-colour diagram so that only those present 
in all jittered frames are used. Only those objects with a photo- 
metric error of 0.207 mag or less in which is our l-cr error in 
K at 19 mag, are used. Fig. |3] shows the colour-colour diagram 
constructed from our JHK observations. The continuous and the 
short dashed green curves show the loci of the colours of main- 
sequence and giant stars adopted from Tokunaga (12000b . The 
long dashed green line shows the loci of Classical T Tauri stars 
(CTTS) from Meyer, Calvet & Hillenbrand (1997 ). The dotted 
red arrows show the reddening vectors up to Ay =50 derived 
from the interstellar extinction law given in Rieke & Lebofsky 
(119851 ). The objects located in between the two main- sequence 
reddening vectors are expected to be reddened main- sequence 
stars. YSOs occupy the region of the colour-colour diagram 



below the lower reddening vector for the main-sequence stars, 
where CTTS, Herbig Ae/Be (HAeBe) stars and Luminous YSOs 
occupy diflferent regimes (Lada & Adams fT992l) . 

Fig. |4] shows the WFCAM ^-band image of the 5.5^x5.5' 
field shown in Fig. [T] All sources in this field exhibiting IR ex- 
cess in the JHK colour-colour diagram (Fig. [3]) are enclosed in 
green circles; sources with IR excess in regions "CI-CIV" are 
enclosed in red boxes. Contours are overlaid on Fig. [4] to show 
the number of objects with IR excess, within a radius of 10'', 
estimated at T'intervals over the 5.5' x5. 5' field. The contours 
show enhancements in the number density of sources with IR 
excess towards regions "CI-CIV" suggesting that these regions 
host embedded clusters. 

The WISE 3.4-22-yum images provide a deeper view of the 
star forming activity in this region. The 12-yum image (Fig. [5]) 
shows a cloud of ~1 5' x30' extending NS. The 22-yum image also 
shows the nebulosity from the whole cloud, albeit with lower 
S/N than in the 12-yum image. The cloud appears cometary; most 
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Fig. 3. The JHK colour-colour diagram of the 13.5 ' x 13.5 ' field surrounding IRAS 07422. The curved lines show the loci of main 
sequence (continuous, green) and giant (short-dashed, green) stars from Tokunaga( 2000 ). The long-dashed green line shows the loci 
of CTTS from Meyer et al. (119971) . The dotted arrows (red) show the reddening vectors for Ay < 50 derived from the interstellar 
extinction law given by Rieke & Lebofsky ( 1985 ). The black dots show objects with photometric error less than 0.207 mag, which 
is our l-cr error for a source of 19 mag in K. Sources which are fainter than 19 mag in K are not shown here. Objects with upper 
limit magnitudes in / are shows with an upward directed arrow; source #1 with H magnitude also is treated as an upper limit is 
shown with upward and rightward directed arrows. 



of the star formation activity seems to be taking place in the 
"head" - in a region of ~6' diameter located towards the southern 
region of the cloud (region encompassed by the white contour 
in Fig. [5]). This region of the cloud exhibits extended nebulos- 
ity in all four bands and hosts all the young clusters discussed 
here. Fig.|6]shows a WISE 3.4, 4.6, 12 //m colour composite im- 
age of a ~10'x 10' field around the southern region of the cloud. 
The nebulosity, which is much stronger at 12yum than at shorter 
WISE wavelengths, is seen red in colour. 

The three MSX sources detected near IRAS 07422 are lo- 
cated in "CI", "CII" and "CIV" (close to the near-IR sources #1, 
#8 and #11 respectively). The three sources detected by AKARI- 
IRC at 1 8 yum agree with the MSX sources within positional er- 
rors. The AKARI-FIS detected emission from "CI" and "CII" 
in all four bands. The FIS source associated with "CI" is the 
brighter one. 



The outflows detected in H2 and the interesting near-IR 
sources detected are described throughout the text. The coor- 
dinates and magnitudes of the prominent sources are given in 
Tabled 

3. ■/. The outflows detected and their possible driving sources 
3.1.1. CI 

"CI" is the most prominent feature in the near-IR images (Fig.[T]). 
This region hosts a cluster of red sources embedded in nebulos- 
ity. An expanded view of a TxTfield containing "CI" is shown 
in Fig.|7] A few point sources of interest are labelled "1-7" on 
the figure. The brightest of these are the two red near-IR sources 
labelled #1 and #2. In the WFCAM images, #1 appears to be 
located behind an extinction lane running NE-SW. Strong nebu- 



5 



Watson R Varricatt: Multiple outflows from the luminous YSO IRAS 07422-2001 




Fig. 4. WFCAM ^-band image of the 5.5'x5.5' field shown in Fig.[lj All objects in this region exhibiting IR excess in the JHK 
colour-colour diagram (Fig. |5]) are enclosed in green circles. The sources in regions "CI-CIV" are enclosed in red boxes. The 
contours in black, yellow, cyan, blue and white respectively enclose regions with 10, 8, 6, 5 and 4 sources with IR excess in a circle 
of lO^'radius, estimated at intervals of T'. 



losity is seen in the vicinity of #1; the nebula is roughly bipolar 
in appearance in / and H bands with the two blobs of diff'use 
emission located NW and SE of the location of #1 and of the 
extinction lane. /, H and K images of the central region of "CI" 
containing #1 and #2 are shown separately as insets on Fig. [7] 

The point source #1 is bright in ^ (10.52 mag) and is much 
fainter in H (14.51 mag). In H band, the emission from the SE 
nebula is stronger than that from the point source detected. In /, 
#1 is not detected and only the nebulosity associated with it is 
seen. Since #1 is closely associated with nebulosity, we derived 
its H and K magnitudes by performing aperture photometry in a 
1.5 ''-diameter aperture. #1 is not detected in /, so we have used 
the / detection limit to derive its / - // colour. The H magnitude 
derived may still be significantly aff'ected by nebulosity. Hence 
we treat the H magnitude also as an upper limit only for deriving 
its colours. Fig. [3] shows that #1 is the reddest object detected in 
this region. #2 is located ~8.85''S of #1. #2 shows large redden- 
ing and excess, but is much less red than #1. There are several 
other sources in "CI", which exhibit large reddening and excess 



in Fig. [3] "4-7" are the prominent ones among them. Prominent 
sources in "CI" with IR excess are shown by "star" symbols in 
Fig. [3] 

The brightest of the three AKARI-IRC sources is located in 
"CI", between #1 and #2; the brighter of the two AKARI-FIS 
sources is in "CI". "CI" also hosts the brightest of the three MSX 
sources, which is detected in all four MSX bands. The location 
of the brightest of the sources detected by IRC, FIS and MSX 
coincides with the IRAS source within the positional errors. 

Fig. [8] shows the 11.2-yum image obtained using UKIRT and 
Michelle. An 1 8'' x 18'' field, which encompasses the objects de- 
tected at 11.2-yum, is shown in the figure. Only three objects 
are detected in the Michelle image; #1 and #2 are the brightest. 
Astrometric corrections are applied to the Michelle image, us- 
ing the coordinates of #1 and #2 from our WFCAM images. The 
third source (labelled "3") is very faint; it is located 3.8 arcsec 
NE of #2. The coordinates of #1 and #2 from the WFCAM im- 
ages and that of #3 measured from the calibrated Michelle im- 
age are listed in Table |2] The 11.2-yum flux obtained from the 
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Table 2. A few of the interesting sources detected 





Coordinates (J2000) 




J 




K 




1 1 2-/ym flux 




(a, 6) 
















1^ 


07-44-27 964 


-20:08 


:30.98 




14.511 (0.002) 


10.516 


(0.020) 


579 (31) 


2^ 


07:44:27.988, 


-20:08 


:39.77 


16.475 (0.003) 


13.317(0.001) 


10.745 


(0.020) 


364 (25) 


3 


07:44:28.251, 


-20:08 


:39.28 










80 (15) 


4b 


07:44:28.196, 


-20:08 


:39.19 




17.933 (0.031) 


14.905 


(0.010) 




5 


07:44:28.552, 


-20:08 


:39.30 




17.463 (0.016) 


15.149 


(0.013) 




6 


07:44:28.123, 


-20:08 


:46.77 


17.227 (0.030) 


15.120 (0.001) 


13.575 


(0.006) 




7 


07:44:27.457, 


-20:08 


:20.92 


16.361 (0.002) 


14.005 (0.004) 


12.119 


(0.002) 




8^ 


07:44:38.814, 


-20:07 


:42.27 


16.920 (0.014) 


13.865 (0.004) 


11.103 


(0.002) 




9b 


07:44:39.035, 


-20:07 


:41.31 


19.360 (0.051) 


17.003 (0.018) 


14.919 


(0.010) 




10 


07:44:35.411, 


-20:07 


:52.88 


19.351 (0.053) 


17.293 (0.003) 


15.125 


(0.002) 




11 


07:44:29.476, 


-20:09 


:57.95 


15.425 (0.015) 


14.271 (0.001) 


12.701 


(0.002) 





^ The values given in parenthesis are the 1-cr errors. 

^ The JHK magnitudes for these sources are aperture photometry in a 1.5 ''-diameter aperture; there may still be contamination from nebulosity. 



Michelle observations are also listed in the table. #1 is the bright- 
est source in the Michelle image at 1 1 .2 jim and #3 is the faintest. 
#3 is detected above 5cr at 1 1.2 micron. We do not see any JHK 
counterpart at exactly the same location as of #3. However, Fig. 
[T] shows a source labelled #4 located only O.IT'W of #3. This 
oflfset is significant, considering the positional accuracy of our 
WFCAM observations. Contours generated from our WFCAM 
^-band image are overlaid on the Michelle image in Fig. [8] 
which show the detection of #1 and #2 at 1 1.2 yum and the oflfset 
of #3 from #4. However, since the Michelle image is astrom- 
etry calibrated with the coordinates of only two sources from 
WFCAM (#1 and #2), we need more observations, especially at 
thermal IR to ascertain if #3 and #4 are the same or diflTerent. #4 
is a highly reddened source. It is detected well in K, marginally 
detected in H and is not detected in /. 

Our continuum-subtracted H2 image (Fig. [2]) reveals several 
line emission features in the vicinity of the cluster. These are la- 
belled MHO 1412-1417 on the figure. The possible directions 
and lengths of some of the outflows implied by these line emis- 
sion knots are shown by dashed arrows on the figure and are 
listed in Table |3] The most prominent of these is MHO 1412, 
which appears to be directed away from #1, with H2 knots ex- 
tending up to a projected length of 27'X0-53 pc) from #1. With 
#1 being the reddest and the brightest source seen in the near- 
and mid-IR, and being the primary contributor to the emission 
seen at longer wavelengths, we identify #1 as the most luminous 
YSO in "CI". From the locations of the features MHO 1412- 
1417, it appears that at least 6 diflferent outflows may be re- 
sponsible for them. Other than #1, we identify #2 and #3 as 
two other YSOs responsible for driving the outflows in this re- 
gion. However it is difl&cult to assign any of the observed H2 line 
emission features MHO 1413-1417 as originating from these 
sources. Some of the other reddened IR-excess sources in "CI" 
may also be driving outflows, some of which are seen in the H2 
line emission. The spectral energy distribution of #1 and a model 
fit to it are discussed in §3.2. 1| 

The extinction lane close to source #1 deserves detailed in- 
vestigation. In / and H bands, it has the appearance of a nearly 
edge-on disk at an angle of -53° east of north (see the inset in 
Fig. [7]). It is also noteworthy that the two H2 knots in MHO 1416 
and the bright knot at the base of MHO 1414 align well at an 
angle of 127° east of north (see Fig. |2]), which is roughly per- 
pendicular to the orientation of the extinction lane. High angular 
resolution observations at longer wavelengths, especially in the 



Table 3. Outflow lengths and directions 



Outflow 


Angle n 


Len| 




ID 


(E ofN) 


arcsec 


parsec 


MHO 1412 


353 


27.2 


0.53 


MHO 1414 


38 


>14.5 


>0.28 


MHO 1415 


41 


>15.5 


>0.30 


MHO 1418 


116 


>13 


>0.25 



mm, are required to understand if the extinction lane is caused 
by a nearly edge-on circumstellar disk. 

3.1.2. CM 

"CII" hosts an embedded cluster -2.7' NE of the IRAS posi- 
tion (Fig.[T]). An expanded view of "CII" is shown in Fig. [9] The 
WFCAM image shows a large number of objects surrounding a 
bright source (labelled #8) detected in the near-IR and exhibiting 
reddening and excess in the JHK colour-colour diagram (Fig. [5]). 
#8 is the brightest near-IR source in "CII". It is less embedded 
than #1 and is detected well in our /, H and K images. There is 
no IRAS detection in "CII". The second brightest of the AKARI- 
IRC and FIS, and MSX sources is detected in this region, but 
they lack the positional accuracy to assign the detection to #8 or 
any other reddened source in the cluster. However, the WISE im- 
ages with very good positional accuracy shows a bright source, 
the position of which agrees well with what is derived from the 
WFCAM images. Hence, #8 is likely to be the most luminous 
YSO in the cluster. In addition to #8, "CII" hosts several other 
reddened sources with IR excess. Prominent sources with ex- 
cess, detected well in our images are shown with blue circles in 
Fig. [3] 

Fig. [2] shows aligned H2 line emission knots labelled 
MHO 1418 and shown by a dotted arrow. The direction of orien- 
tation of the H2 knots imply an outflow emanating from "CII". 
It is not clear if #8 or any of the fainter reddened sources with 
excess nearby is responsible for driving MHO 1418. The limited 
H2 knots in MHO 1418 do not trace directly to #8, but they trace 
back to a source #9 located 3.25''NE of #8. MHO 1418 has a 
length of W\03l pc) if it is from #8 and 13'' (0.25 pc) if it is 
from #9. Fig. |3] shows that #9 also exhibits large reddening and 
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Fig. 5. The WISE 12-yum image over a field of 18'x33' encom- 
passing the cloud in which IRAS 07422 is located. The contours 
in black, white and grey show the boundaries where the emis- 
sion from the cloud falls to 15%, 40% and 85% respectively of 
the peak flux in the bar. Most of the star formation happening 
within this cloud is in the region within the white contours. 



excess. #8 is the dominant source in the WISE bands. In ^3.2.2[ 
we model the spectral energy distribution of #8. 

The WISE 12-yum image shows a nebulous emission feature 
at Qf=07:44:40.034, ^=-20:08:15.91, labelled "Blob" in Figs.|5] 
and[6] It is off'set ~39''SE of #8. The WISE 3.4 and 4.6-yum im- 
ages show a point source very close to this position. It is not clear 
if they are the same as "Blob". The 22-yum image has much poor 
spatial resolution. It doesn't resolve a source at the location of 
"Blob", but only shows a slight elongation of the nebulosity as- 
sociated with #8 in the in the direction of the "Blob". Therefore it 
remains to be investigated at better spatial resolution and higher 
sensitivity at longer wavelengths to learn if "Blob" hosts a source 
much younger than #8. 



^(CII) 



Blob 5- 



• 



11 • 

m 



Fig. 6. The WISE colour-composite image (3.4 yum-blue, 4.6 jim- 
green, 12yum-red) in a 10' x 10' field (the region enclosed in box 
in Fig. [5]) including the embedded clusters and the bright region 
of the cloud containing them. 




Fig. 7. WFCAM JHH2 colour-composite image (/-blue, H- 
green, H2-rQd) in a VxV field containing "CI". The top-left inset 
shows /, H and K images of the central 1 2'' x 18'' field (enclosed 
in the dashed box) at a higher contrast, with circles showing the 
locations of the two bright sources derived from the ^-band im- 
age; a dotted yellow line is drawn on the /-band image along 
the extinciton lane to guide the eyes. The top-right inset shows 
a JHK colour composite of the same region (/-blue, i/-green, 
^-red), revealing the bright point sources embedded in the neb- 
ulosity. 
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Fig. 8. Michelle image of IRAS 07422 at 11.2yum covering a 
field of 18''xl8''. Contours from our WFCAM /^-band image 
are overlaid. The "+" shows the location of the MSX source. 



are several other reddened sources in "CIII". In H2 line, we de- 
tect three line emission knots, which are labelled MHO 1419 in 
Fig.|2] It is not clear if all these features are from YSOs in "CIII". 
One of these knots is is located very close to #10; it may be as- 
sociated with #10. With the WISE source coinciding with #10, 
we identify it as the dominant YSO in "CIII". 

"CIV" is located towards the south of "CI"; it also hosts 
several reddened sources as seen in Figs. [T] and [3] Sources in 
"CIV" with excess are shown with cyan boxes in Fig. |3] The 
source labelled #11 is the prominent among them and is de- 
tected well in /, H and K. Other than "CI" and "CII", "CIV" 
is the only region in Fig. [T] detected by MSX and AKARI-IRC. 
Most of the reddened sources in "CIV" except #1 1 appear to be 
distributed in a filamentary pattern oriented NE-SW. #1 1 is lo- 
cated somewhat south of those reddened sources. At the spatial 
resolution of WISE, more than one source may be contributing 
to the WISE detection at the location of #11. However the lo- 
cation of the WISE source shows that #1 1 is the dominant con- 
tributor to the flux measured in the WISE bands. The faintest of 
the three AKARI-IRC sources coincides with the near-IR source 
#11. This object was not detected by AKARI-FIS in the far-IR. 
This is also the faintest of the three sources detected by MSX 
in this region. "CIV" was detected by MSX in band A only. We 
detect an H2 knot MHO 1420, NE of "CIV". It is not clear if it 
is due to outflow from any of the sources in "CIV" or "CI" dis- 
cussed here. The WISE images exhibit a bar-like feature in the 
nebulosity below "CIV", which is labelled in Figs. [5] and [6] 



■ % 



Fig. 9. WFCAM JHH2 colour-composite image (/-blue, H- 
green, if2-red) in a 1.5^x1.25' field containing "CII" and "CIII". 



3.1.3. CIII and CIV 

Regions "CIII" and "CIV" show concentrations of red sources. 
These appear elongated and are much less prominent than "CI" 
and "CII". Fig. [4] shows that there is enhancement in the number 
density of sources with IR excess towards "CIII" and "CIV". 
They are shown enclosed within ovals in Fig. [T] Deep thermal 
imaging at high angular resolution is required to understand if 
these two regions host embedded clusters. 

"CIII" is seen elongated in the NE-SW direction and is lo- 
cated close to, and towards the west of, "CII". Several sources 
in "CIII" show IR excess, which are shown with green triangles 
in Fig.|3] There is a very red source, #10, located somewhat near 
the centre of "CIII", and showing large reddening and excess in 
Fig. [3] It is detected well in H and K bands and is marginally 
detected in /. #10 is detected well in all four WISE bands; in the 
WISE 22-yum image, it is one of the few objects detected. There 



3.2. The spectral energy distribution 

Out of the four embedded clusters identified in our data, "CI" 
and "CII" are the prominent ones. Near-IR sources #1 and #8 ap- 
pear to be the most luminous YSOs in these clusters, as inferred 
from their possible association with longer wavelength data. The 
SEDs of these sources are modelled in ^3.2.1| and ^3.2.2| 

The WFCAM, Michelle, MSX, WISE, AKARI and IRAS 
data were used to construct the SED of the sources. Colour cor- 
rections were applied to AKARI and IRAS data using the cor- 
rection factors given in their respective point source catalogues 
(Kataza et al. 120101 Yamamura et al. 12010 1, Beichman et al. 
fT988l) . 

We used the SED fitting tool of Robitaille et al. (120071) avail- 
able on-line to model the SEDs. They use a grid of 2D radia- 
tive transfer models presented in Robitaille et al. ( 2006), devel- 
oped by Whitney et al. (2003a, 2003b, etc.). The grid consists of 
20000 YSO models with SEDs covering stellar masses of 0.1- 
50 M© and evolutionary stages from the early infall stage to the 
late disk-only stage, each at 10 diflferent viewing angles. The re- 
sults of the SED modelling are presented in Table [4] Note that 
the model assumes a single object whereas multiplicity will be 
an issue in our regions. Hence, the actual errors in some of the 
parameters could be larger than the values given in Table [4] due 
to possible contributions to the observed fluxes of the dominant 
YSOs from neighbouring sources. 

3.2.1. Source #1 

At the poor spatial resolution of IRAS, AKARI-FIS and MSX, 
the source detected in "CI" will contain the emission from the 
whole cluster. However, #1 is resolved from the other sources in 
"CI" in the WISE 3.4 and 4.6-yum data. The WISE images and 
the well-resolved Michelle image, both with good positional ac- 
curacy, show that #1 is the dominant contributor at longer wave- 
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lengths (see Fig. [8]). #1 has a steeply rising SED. In K band, the 
combined emission from #2 - #7, the extended nebulosity and 
the rest of the objects in the vicinity dominates over the emission 
from #1; the relative contribution from #1 being only ~21%. At 
11. 2 yum, #1 contributes 56.6% of the total flux of #1, #2 and #3 
combined. 

For //, K and 11.2yum, we used the flux estimated from the 
WFCAM and Michelle data (see ^3.1.1 ) for constructing the 
SED. The WISE photometry of #1 is quite free of emission from 
the neighbours, especially in Wl and W2 bands, where #1 is re- 
solved from #2 and the rest of the cluster members. However, 
there appears to be a small contribution from the neighbours in 
the photometry listed in the catalogue. Hence, for WISE Wl and 
W2 bands, we re-determined the photometry of source #1 by 
performing relative photometry using the magnitudes of neigh- 
bouring isolated point sources in the WISE images and using 
a IC -diameter aperture. In WISE band W3, #1 is not resolved 
from #2 and the rest of the members of the cluster due to poor 
spatial resolution. We therefore estimated the photometry of the 
whole cluster using a large aperture and applied a flux correction 
factor of 0.57, which is the relative contribution of #1 towards 
the flux of #1, #2 and #3 combined, derived from the Michelle 
data at 11.2yum given in Table |2] The IRAS 12yum flux has also 
been scaled by the same factor. No correction has been applied 
for the WISE 23-yum photometry assuming that emission at long 
wavelengths is dominated by the emission from #1. 

The flux in MSX band A includes emission from the rest 
of the cluster. MSX flux in band A is corrected by the relative 
contribution from source #1 towards the observed flux, which 
is the emission from the whole cluster. This correction factor is 
derived by interpolating for the value at 8.28 yum between the 
relative contribution of #1 in the WISE 4.6-yum data and in the 
Michelle 11.2-yum data. No corrections are applied to the MSX 
band E data. MSX fluxes of #1 in bands C and D are not used in 
our fit since these are of poor quality. 

The IRAS lOOyum flux of this region is of moderate quality 
only and will be contaminated by additional emission from the 
cluster at the low spatial resolution of IRAS. It is also aff'ected 
by IR cirrus. So it is not included in the model fit. 

The left panel in Fig. [TO] shows the model fit to the SED of 
source #1. Column 2 of Table |4] shows the results of the SED 
modelling for source #1. In addition to the best fit model, there 
were only two additional models with (x^ -xlestfu^ ^^^^ point 
< 3. Out of those, one diff'ers from the best fit model only in fore- 
ground and circumstellar Ay and the angle of inclination of the 
disk axis. Hence, the errors given against the parameter values 
in Table |4] are the diff'erences in the parameter values of the best 
fit model and of the additional fit. 

We derived a mass, luminosity and temperature of 10.7 M©, 
3.9x10^ L© and 7840 K respectively from the fit. #1 is actively 
accreting mass and it is likely to be a Class-I protostar. A 
foreground visual extinction Ay=16mag is obtained from the 
model fit, whereas Fig. |3] implies a much larger Ay for the 
IR source. This suggests that #1 has a large amount of cir- 
cumstellar matter, probably as a disk. This is implied by an 
Av(circumstellar)=25.4mag given by the fit. High angular res- 
olution observations in the far-IR to sub-mm wavelengths are 
required to further constrain the model fit. 



3.2.2. Source #8 

WFCAM JHK, AKARI-IRC and FIS, WISE and MSX band-A 
data are used for fitting the SED of #8. MSX C, D and E-band 



data are not included in the final fit since they are of poor quality; 
(nevertheless, they agree well with the current fit and including 
them in the fit did not change the fit results). WFCAM JHK 
photometry are obtained in a 1.5 ''-diameter aperture to reduce 
the contamination from the neighbours and the nebulosity in the 
cluster to a minimum. 



The right panel in Fig. 10 shows the model fit to the SED 
of #8. Table |4] shows the parameters of the best fit model. Upon 
the initial attempts to fit the data, it was noticed that AKARI- 
FIS data beyond 65 yum exhibit excess with respect to the model 
fit, which systematically increases towards longer wavelengths. 
Hence the AKARI 90, 140 and 160-yum data were not included in 
the fit and were treated as upper limits with low probability only. 
These data points are are shown as downward directed triangles 
in the figure. The data points shortward of 65 yum fit the model 
well; a stellar mass, temperature, luminosity and age of 6.7 M©, 
19950 K, 1.4x10^ L© and 6.0x10^ yrs respectively are derived. 
#8 appears to be older and much less embedded than #1. The 
model fit yields a foreground extinction similar to that for #1, 
but a much lower circumstellar extinction. The high temperature 
of #8 suggests that it might have formed an Hii region around 
it. However, note that even after formation of an Hii region, the 
young star can still grow through accretion of ionized gas (eg. 
Keto|2007l). 

Nevertheless, it should be noted that the temparature of #8 
derived by the model fit (19950 K) is too high even for a main 
sequence star of its derived mass of 6.7 M©. This can happen 
if there is extra contribution to the flux observed at short wave- 
lengths which may make the object appear hotter. It should also 
be noted that a circumstellar extinction of Ay=0.64mag is ab- 
normally low for a young source exhibiting a large IR excess as 
seen in Fig. [3] It is very much possible that either an angle of 
49.5° estimated by the SED modelling is wrong or more than 
one source is contributing to the near- and mid-IR fluxes that we 
use for model fit. High angular resolution observations are re- 
quired for #8 to understand the real nature of the dominant YSO 
here. 

The deviations of the AKARI-FIS data from the model fit, 
which increase towards longer wavelengths, suggest an extra 
contribution in the beam at longer wavelengths from a source 
colder than #8. The feature labelled "Blob" (Figs. [5] and |6] and 
^3.1.2 ) needs to be studied at better spatial resolution and higher 
sensitivity at longer wavelengths to see if the it hosts a source 
much younger than #8 and if that is the one contributing to the 
emission in the FIS bands, which causes the FIS fluxes to deviate 
systematically from the model fit. Some other near-IR sources 
like #9, also need to be understood in detail. 



3.3. The general picture massive star formation in 
IRAS 07422 as seen in the infrared 

At a distance of 4 kpc, the cloud detected in the WISE images 
is larger than -17.5 pcx35 pc. Most of the star forming activ- 
ity seems to be happening in a region of diameter <11.5pc in 
the southern region of the cloud. The four clusters/cluster candi- 
dates "CI-CIV" detected by us are in this region. We see H2 line 
emission features which are likely to be caused by jets from "CI- 
CIII". The brightest object in "CIV" also has infrared colours 
typical of YSOs. At least in clusters "CI-CIII", the brightest ob- 
ject in near- and mid-IR seems to be located near the centres of 
the clusters. 

Several theoretical and observational studies of massive star 
forming regions suggest that massive star formation is intimately 
linked to the hierarchical formation of massive stellar clus- 
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Fig. 10. The left panel shows the SED of #1 . The filled circles show the data from WFCAM at H and K, MSX at 8.28 and 21 .34 yum, 
Michelle at 11.2yum, WISE at 3.6-22 yum, IRAS at 12-60 yum, and AKARI at 18-140 yum wavelengths. The right panel shows the 
SED of #8. The filled circles show the data from WFCAM at /, H and MSX at 8.28 yum, Michelle at 1 1.2 yum, WISE 3.6-22 yum 
and AKARI at 18 and 65 yum; the downward directed triangles show the AKARI 90-160yum-band data, which are treated as upper 
limits only. The continuous line shows the best fit model and the grey lines show subsequent good fits for (x^ - xlestfit^ ^^^^ 
point < 3. The dashed line corresponds to the stellar photosphere for the central source of the best fitting model, as it would look in 
the absence of circumstellar dust (but including interstellar extinction). 



Table 4. Results from SED fitting 



Parameter 


Best fit values^ 




Source #1 


Source #8 


Stellar mass (M©) 


10.7 (+0, -0.2) 


6.7 


Stellar age (yr) 


3.0 (+0, -0.2)xl04 


6.0x10^ 


Stellar radius (R©) 


33.6 (+4.5, -0) 


3.1 


Stellar temperature (K) 


7840 (+0, -625) 


19950 


Disk mass (M©) 


1.3 (+0, -0.3)xlO-2 


3.6x10-3 


Disk accretion rate (Mgyr"^) 


2.9 (+0, -2.42)xl0-6 


7.6x10-^ 


Disk/envelope inner radius (AU) 


8.4 (+2.2, -0) 


12.7 


Disk outer radius (AU) 


16.1 (+1.2,-0) 


134 


Envelope mass (M©) 


6.1 (+0.5, -0)xl02 


2.3 


Envelope accretion rate (M©yr~^) 


9.3 (+l.l,-0)xl0-4 


7.5x10-^ 


Envelope outer radius (AU) 


1.0x10^ 


4.3x10^ 


Total Luminosity (L©) 


3.9 (+0.3, -0)xl03 


1.4x10^ 


Angle of inclination of the disk axis^ (°) 


31.8 (+18.5) 


49.5 (+14.6, -8.1) 


Ay (Foreground) 


16.1 (+4.7,-0.8) 


17.1 (+0, -0.5) 


Ay (Circumstellar) 


25.4 (+0, -8.7) 


0.6 (+0, -0.1) 



^ The values given in parenthesis are the differences in the parameters of the additional models with (x^ ~xlestfi) P^^ ^^^^ point < 3, from those 
of the best fit model. No errors are given when the parameters of these models do not diff'er; nevertheless, the uncertainties could be higher as 
described in |3.2| A distance of 4.01 kpc is used. 

^ The angle of inclination of the disk axis is with respect to the line of sight. 



ters, where massive stars form near the centres of massive sub- 
clusters and grow by competitive accretion aided by the cluster 
potential. In this scenario, several subclusters form in a massive 
clump, massive stars near the central regions of these subclus- 
ters grow through competitive accretion and finally the subclus- 
ters m erge (eg. Testi et al. 2000, Clarke, Bonnell & Hillenbrand 
12000', Bonnell, Bate & Vine 2003 , Bonnell & Bate 2006 , Bonnell 
& Smith 2011 ). Our observations reported here support this sce- 
nario. (However, note that precise radial velocity studies to un- 
derstand the dymamics of the gas associated with this region 
are required for a confirmation). At least in "CI" and "CII", the 



brightest objects in the infrared are seen close to the centres of 
the clusters hosting them. In "CIII" also, the brightest object 
seems to be located close to the centre of the elongated pat- 
tern of distribution of the reddened objects. WISE observations, 
though deep, lack the spatial resolution to resolve the members 
of the clusters at these large distances. More sensitive ^-band 
and thermal IR imaging with high angular resolution are re- 
quired to understand the census and distribution of the young 
stars in these clusters. High angular resolution observations at 
sub-mm and mm wavelengths are also required to understand if 
there are younger members in these clusters than those we see at 
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near-IR and to get a more detailed picture of the star and cluster 
forming activity in the cloud. 

4. Conclusions 

1. At least two (and possibly four) young subclusters, embed- 
ded in a cloud are discovered in IRAS 07422. Narrow-band 
H2 1-0 S(l) line images reveal eight or more outflows from 
them. 

2. The most prominent and youngest cluster in this region is 
"CI", associated with the IRAS source and the brightest of 
the MSX and AKARI sources. Multiple star formation is 
happening in "CI"; H2 line emission features observed in its 
vicinity are suggestive of at least six outflows produced by 
YSOs in "CI". #1 is the most luminous object in this region. 
We derive a mass and luminosity of 10.7 M© and 3.9x10^ Lo 
respectively by modelling its SED. This object appears to be 
very young and and is probably a Class-I protostar. High an- 
gular resolution observations at 11.2^m reveal three of the 
embedded sources in "CI". 

3. Source #1 in "CI" is located close to an extinction lane, in- 
clined at -53° east of north. A set of three H2 knots are seen 
aligned in a direction roughly perpendicular to the extinction 
lane. High angular resolution at mm wavelengths are reqired 
to understand if the extinction lane is a circumstellar disk 
seen nearly edge on. 

4. #8 is the brightest object in the near-IR, located in the sec- 
ond most prominent cluster in this region, "CII". The second 
brightest object in the AKARI, WISE and MSX appears to 
be associated with this region, but the spatial resolution of 
any of these surveys is not suflficient to identify between #8 
and any of its very close neighbours. There is an outflow de- 
tected from this cluster. It is not clear if #8 or any of it very 
close neighbours is responsible for the outflow. Modelling 
the SED of #8 yields a mass, luminosity and temperature of 
6.7 Mo, 1.4x10^ Lo and 19950 K respectively. This temper- 
ature is high even for a main sequence star of 6.7 M© sug- 
gesting possible extra contributions to the flux at short wave- 
lengths. The SED fit shows the four AKARI-FIS fluxes to be 
deviating systematically from the fit, with the residuals in- 
creasing towards longer wavelengths. These are suggestive 
of contribution from a colder component to the flux mea- 
sured by AKARI. Cluster candidates "CIII" and "CIV" also 
host YSOs; "CIII" appears to host one or more outflows. 

5. Multiple star formation is taking place inside more than one 
sub-cluster within a cloud in this region. This region appears 
to be a case of hierarchical mode of cluster formation where 
massive star formation takes place in multiple sub-clusters in 
a cloud and the massive stars at the centres of the sub-clusters 
grow through competitive accretion. 

6. Observations at IR through mm wavelengths at better angu- 
lar resolution are required for a more detailed understanding 
of the star formation activity in this cloud. 
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